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FOREWORD 

This  report  was  prepared  under  Work  Order  No.  086432  of  Contract  No. 
N0001 4-78-204  in  support  of  the  U.S.  Naval  Postgraduate  School  research 
project  supported  by  the  Naval  Air  Systems  Command  (AIR  370)  and  the  Naval 
Material  Command  (EO/MET).  The  data  were  obtained  by  the  Environmental 
Physics  Group  at  NPS  under  the  direction  of  Professors  K.  L.  Davidson  and 
G.  E.  Schacher. 
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ABSTRACT 


''^Aerosol  and  Meteorological  data  from  cruises  in  the  East  Pacific 
(CEWC9M-78)  and  the  North  Atlantic  ( JASIN)  were  used  to  evaluate  the 
predictive  performance  of  the  Wells-Munn-Katz  aerosol  model.  Under 
given  meteorological  conditions  (wind  speed  and  relative  humidity), 
the  model  gave  very  good  predictions  of  the  average  aerosol  extinction 
coefficient  at  10.6jCm  wavelength  (the  standard  deviation  was  about  a 
factor  of  three).  The  model  did  not  perform  well  when  asked  to  predict 
specific  aerosol  extinction  values.  The  model  continental  coefficient, 

B  =  1.7,  was  found  to  be  too  large-T^the  average  open  ocean  value  should 
be  B  =  0.24.  The  continental'  coefficient  was  found  to  be  important  for 
predicting  IR  band  extinction  from  visible  band  extinction  estimates. 
Further  analysis  demonstrated  that  much  of  the  variance  of  continental 
aerosol  spectral  density  was  due  to  air-mass  history  and  that  much  of  the 
variance  of  sea  salt  aerosol  spectral  density  was  due  to  changes  in  the 
marine  layer  mixing  height. 
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A.  INTRODUCTION 


Aerosol  size  spectrum  data  from  two  oceanic  field  experiments  have 
been  used  to  evaluate  the  Wells-Munn-Katz  (WMK)  marine  aerosol  Model  B 
(Noonkester,  1980).  The  field  experiments  were  the  Cooperative  Experi¬ 
ment  for  West  Coast  Oceanography  and  Meteorology  (CEWC0M-78)  off  the 
California  coast  in  May,  1978  and  the  Joint  Air-Sea  Interaction  (JASIN) 
experiment  between  Scotland  and  Iceland  in  August  and  September,  1979 
(see  Figure  1).  The  purpose  of  this  report  is  to  evaluate  the  accuracy 
of  WMK  predictions  of  aerosol  extinction  coefficient,  a,  and  to  produce 
modifications,  if  possible,  that  improve  the  model. 

B.  MEASUREMENTS 

1.  Aerosols 

The  aerosol  size  spectra  were  compiled  using  Particle  Measure¬ 
ments  Systems,  Inc.,  (PMS)  optical  spectrometers.  The  smaller  sizes 
(0.075  to  1.5um  radius)  were  counted  with  a  PMS  model  ASAS-100  (active 
scattering)  probe;  the  larger  sizes  (0.2  to  14um  radius)  were  counted  with 
a  PMS  Model  CSAS-300  KV  (classical  scattering)  probe.  The  spectrometer 
probes  were  controlled  by  a  PMS  Model  DAS-32  data  acquisition  system  with 
the  data  recorded  on  digital  magnetic  tape.  Further  details  on  the 
aerosol  system  can  be  found  in  Schacher  et  al  (1980). 

2.  Meteorology 

The  meteorological  parameters  of  direct  relevance  to  the  WMK 

model  are  relative  humidity  (H)  and  true  wind  speed  (u).  In  addition,  we 

obtained  measurements  of  air  temperature  (T),  sea  surface  temperature 

(T  ),  relative  wind  speed  (UJ  and  direction  (0),  ship  speed  (u  )>  and 
s  r  5 

mixed  layer  height  (h).  A  complete  description  of  the  equipment  is  given 
in  Houlihan  et  al  (1978).  A  summary  is  provided  in  Table  I.  Note  that 
CEWCOM-78  had  more  accurate  T  and  H  measurement,  and  tv/o  complete  levels 
of  meteorological  instrumentation.  The  ship's  speed  (Us)  was  taken  from 
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a  shipboard  instrument  during  JASIN  but  was  simply  "eyeball"  estimated 
by  the  captain  for  CEWCOM-78.  Thus,  the  true  wind  data  are  more  accurate 
for  JASIN. 

C.  ANALYSIS 

1 .  Aerosol  Spectra 

The  basic  aerosol  data  is  obtained  as  N(r)  =  dN/dr  in  90  size 
channels  (6  ranges  of  15  channels  each)  from  0.075  to  14pm  radius.  The 
N(r)  spectrum  is  calculated  for  successive  20  minute  averaging  periods. 

The  first  few  size  bins  of  each  range  are  thrown  out  becuase  they  are 
influenced  by  noise  and  bias  errors  in  the  photodetector.  Other  bins 
in  the  0.5  to  2pm  radius  area  are  thrown  out  because  they  fall  in  the 
"ambiguity  zones"  associated  with  non-single  valued  nature  of  the  light 
scattering  at  the  0.63pm  laser  wavelength  of  the  spectrometers. 

The  N(r)  spectrum  is  fit  in  LGT  (N(r)),  LGT(r)  space  with  a 
seventh  order  polynomial  for  0.1  <r  <7ym  and  a  first  order  polynomial 
(linear)  for  r>  7ym  .  The  linear  fit  (in  log-log  space)  is  used  for  the 
larger  sizes  to  avoid  the  instabilities  associated  with  extrapolating  high 
order  polynomials  beyond  regions  of  poor  statistics.  The  extinction  is 
calculated  using 

30 

a  =  irr2N(r)  E(n,A)dr  (1) 

J0.03 

where  E(n,A)  i,s  the  Mie  scattering  efficiency  at  wavelength  A  and  refractive 
index  n.  This  entire  process  is  described  in  great  detail  in  Schacher  et 
al  (1980). 

2.  Model  Calculations 

We  have  used  the  WMK  model  B  as  aiven  by  Noonkester  (1980).  Since 
several  WMK  models  are  available,  we  will  provide  exact  details  here.  The 
size  spectrum  is  given  by 
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N(r)  =  B(r/a)"4  +  1 .62(C]  +  C2V6)(eV 

where  r  =  radius, 

v  =  0.5  for  u  £  4  m/s 
v  =  u-  3.5  for  u  >  4  m/s 
F  =  1  +  ( V/60) 3 
T  =  0.384  -  0.00293  V1 ,25 
Z  =  height  above  sea  surface,  m 
h0  =  scale  height,  m  (800m  for  Z  <  1000m) 
B  =  1.7 

a  =  0.81  exp  0.066S/O  .058  -  S) 

S  =  H/100 

The  other  constants  are  given  below 


-=£  -  8.5(r/a)r 


m/s 


'1 


v  <7 
v  >7 


350 
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1000 

6900 


1.15 

0.29 
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r 

Fa 


(2) 


We  have  used  the  model  in  this  form  for  the  extinction  calculating 

although  it  is  theoretically  wrong  by  a  factor  a-1  (see  Wells  et  al ,  1977, 

and  Hughes,  1980).  The  first  item  in  Eq.  2  is  the  continental  aerosol 

-4 

component  and  has  an  r  dependence  generally  referred  to  as  a  "Junge" 
distribution.  The  second  term,  which  is  wind  speed  dependent,  is  the 
locally  generated  sea  salt  component. 

3.  Volume  Representation 

Since  the  size  dependence  of  aerosol  spectral  density,  N(r),  is 

_3 

typically  on  the  order  of  r  ,  the  volume  spectrum  is  often  more  convenient 
to  deal  with. 


V(r)  =  dV  =  4  ur5  N(r) 
dr  3 

Thus,  we  can  partition  the  marine  aerosol  volume  distribution  into  con¬ 
tinental  (VJ  and  sea  salt  (V,-)  components 
c  ^ 


(3) 
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V ( r)  =  Vc(r)  +  V$(r)  (4) 

where  we  assume 

Vc(r)  =  A'  a3/r  =  A/r  (5) 

3 

The  quantity  A1  =  A/a  =  Ajt  B  is  called  the  normalized  Junge  coefficient 
(A  is  the  unnormalized  Juflae  coefficient). 

4.  Editing 

There  are  950  JASi.i  and  450  rFWC0M-78  aerosol  spectra  (20  minute 
averages).  The  criteria  to  exclude  data  from  the  analysis  are  listed 
below 

1)  Relative  wind  direction  not  within  30°  of  aerosol  probe  intake 

axes. 

2)  Humidity  and  true  wind  speed  data  unavailable. 

3)  Fog  events  (a  >  1  km~^) 

4)  Continental  influence:  Certain  time  periods  in  CEWCOM-78 

were  heavily  contaminated  with  continental  aerosols  due  to  proximity  to  land 
or  off  shore  wind  conditions. 

5)  Insufficient  wind  speed  data  over  previous  twelve  hours: 

The  true  wind  speed  used  in  most  of  the  WMK  comparisons  was  an  average 

of  over  the  previous  twelve  hours.  Previous  investigations  (Katz,  1980) 
have  shown  that  the  long  time  periods  required  for  the  aerosols  to  reach 
equilibrium  with  the  local  surface  generation  precludes  the  use  of  a  simple 
20  minute  average  wind  speed.  We  excluded  data  if  there  was  more  than  a 
three  hour  gap  in  the  wind  speed  record  over  the  previous  twelve  hours. 

This  eliminated  about  30%  of  the  JASIN  data  and  about  90%  of  the  CEWCOM-78 
data. 
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D.  RESULTS 

The  results  presented  here  will  tend  to  emphasize  the  JASIN  data. 
Because  CEWCOM-78  was  so  close  to  the  continent,  the  variations  in  back¬ 
ground  (continental)  aerosols  was  considerable.  This  problem  was  ac¬ 
centuated  by  a  cruise  plan  that  brought  the  ship  into  the  Los  Angeles 
basin  on  several  occasions.  The  periodic  downwind  cruise  track  of  the 
ACANIA  greatly  reduced  the  amount  of  data  available  for  the  12  hour  wind 
speed  averages. 

One  further  note  on  terminology  is  in  order.  For  the  sake  of 
simplicity,  aerosol  extinction  coefficients  calculated  from  the  aerosol 
spectral  density  data  (Section  C-l)  shall  be  referred  to  as  the  "measured" 
extinction.  Those  who  feel  that  this  is  unbearably  arrogant  are  referred 
to  a  recent  report  (Fairall  et  al ,  1980)  which  documents  the  excellent 
correlation  of  optical  measurements  of  extinction  with  NPS  aerosol 
spectrum  extinction  calculations. 

1.  Model  Extinction  Comparison 

The  basic  measurements  of  extinction  at  three  relevant  wavelengths 
are  presented  in  Figure  2  to  illustrate  the  strong  wind  speed  dependence. 

In  Figure  3  the  data  are  compared  with  the  WMK  predictions  for  A=  10.6  pm 
at  RH  =  87%.  Because  of  the  scarcity  of  CEWCOM-78  data,  we  also  included 
the  20  minute  wind  speed  average  results.  The  JASIN  data  12  hour  averages 
were  a  much  better  fit  to  the  WMK  curve  than  the  20  minute  averages  and, 
in  fact,  were  very  similar  to  the  CEWCOM-78  results.  The  ratio  of  the 
visible  extinction  to  the  IR  extinction  (a  useful  parameter  because  it 
is  essentially  independent  of  relative  humidity)  is  shown  in  Figure  4 
for  10.6  and  3.75pm  wavelengths. 

This  comparison  (Figure  3)  asks  "given  a  wind  speed  and  relative 
humidity,  how  well  does  the  model  predict  the  observed  extinction"? 

Clearly,  the  WMK  model  predicts  very  well  on  the  average  with  rather 
large  standard  deviation  (about  half  an  order  of  magnitude,  ora  factor 
of  3).  A  considerably  more  stringent  comparison  asks  "given  and  observed 
extinction,  how  well  does  the  model  predict  this  extinction"?  This 
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Figure  3a.  Measured  Extinction  Coefficient  (Circles)  and  WMK  Model 

Predictions  at  X  —  10.6ym  and  RH  =  87%  Versus  Wind  Speed. 
J AS i N  (12  Hour  Average  u) , 
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Figure  3b.  Measured  Extinction  Coefficient  (Circles)  and  WMK 
Model  Predictions  at  X  =  10. 6pm  and  RH  =  87 %  Versus 
Wind  Speed.  CEWCOM-78  (Circles  -  12  Hour  Average  u, 
X's  -  20  Minute  Average  u). 
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comparison  (Figure  5)  is  considerably  less  favorable,  particularly  in 
the  extremes  of  very  clear  and  very  unclear  conditions.  A  certain  amount 
of  disagreement  of  this  type  is  expected  when  comparing  quantities  that 
are  subject  to  experimental  error.  However,  we  show  in  Figure  6  that  a 
considerable  part  of  the  discrepancy  is  due  to  over  estimation  of  the 
continental  aerosol  component  (WMK  assume  B  =  1.7  or  A'  =  7.1).  Further¬ 
more,  the  high  visibility  conditions  are  characterized  by  much  smaller 
continental  aerosol  densities  (A1  =  0.2)  while  the  low  visibility  condi¬ 
tions  have  typically  greater  continental  aerosol  densities  (A1  =  1.7). 

The  disagreement  at  high  extinction  values  is  due  to  the  inaccuracy  of  the 
relative  humidity  measurement  which  is  critical  under  near  fog  and  heavy 
haze  conditions  where  RH  =100%. 

2.  Model  Aerosol  Spectrum  Comparison 

Since  the  extinction  results  of  the  previous  section  revealed 
some  problems  with  the  WMK  model,  we  decided  to  use  the  aerosol  spectrum 
volume  density,  V(r),  to  examine  the  model  more  fundamentally.  The  total 
conglomerate  volume  spectra  for  JASIN  are  shown  for  different  12  hour 
wind  speed  categories  in  Fig.  7.  This  graph  clearly  indicates  that  V ( r ) 
is  virtually  independent  of  wind  speed  at  r  =  O.lym.  We  conclude  that 
surface  generated  sea  salt  aerosols  are  insignificant  at  r  =  0.1yn  and  that  these 
particles  are  of  non-local  (continental  background)  origin.  Therefore, 
by  keying  on  the  V(r)  spectrum  for  r<0.3ym,  we  can  simply  calculate  the 
continental  coefficient  from  the  observed  spectrum  using  Eq.  5  (A  =  rV(r)). 

This  method  was  used  to  obtain  the  values  of  A'  given  in  Fig.  6.  Further 
proof  that  this  gives  a  realistic  index  of  continental  influence  is  given  in 
Fig  8  where  A'  (recall  that  A'  is  A  corrected  for  relative  humidity)  is 
compared  to  atmospheric  radon  activity  for  CEWCOM-78  (Larson  et  al ,  1979). 

A  typical  value  of  A'  was  1.0  for  JASIN  with  a  range  from  .1  to 
5.  For  CEWCOM-78  (excluding  several  periods  directly  off  the  coast),  A1 
was  typically  2.0  with  a  range  from  1  to  10.  If  we  simply  change  the 
continental  coefficient  from  A'  =  7.1  (  B=  1.7)  to  A'  =  1  (B  =  0.24), 
the  WMK  model  predictions  are  improved  in  the  clearer  conditions:  (X's  in 
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as  a  Function  of  Extinction  Actually  Observed 
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Fig.  5).  If  we  go  one  step  further  and  use  the  visible  wavelength  ex¬ 
tinction  with  Fig.  6  to  estimate  the  value  of  A',  then  we  see  virtually 
no  improvement  (triangles  in  Fig.  5).  At  this  point,  the  discrepancy  is 
due  to  the  sea  salt  component. 

E.  DISCUSSION 

In  the  previous  section  we  found  that,  if  given  a  wind  speed  and 
humidity,  the  WMK  was  a  good  predictor  of  the  average  extinction  observed, 
although  30%  of  the  observations  disagreed  by  more  than  a  factor  of  3. 

On  average,  the  predictions  based  on  12  hour  averages  of  wind  speed  were 
better  than  those  based  on  20  minute  averages.  The  model  was  less  success¬ 
ful  when  compared  point  by  point  with  the  actual  observed  extinction.  In 
this  case  the  20  minute  averages  were  slightly  better.  Some  of  the  disa¬ 
greement  at  high  visibility  was  corrected  by  using  a  more  realistic  con¬ 
tinental  aerosol  coefficient,  the  remainder  is  due  to  considerable  over¬ 
estimation  of  the  sea  salt  component  in  roughly  20%  of  the  samples. 

It  is  important  to  realize  that  the  WMK  model  is  an  average  continental, 
equilibrium  surface  generation  model .  Deviations  of  the  continental  aerosol 
were  shown  in  Figs.  6  and  8.  The  sea  salt  component  is  also  subject  to 
considerable  deviations  from  equilibrium  (this  is  one  of  the  reasons  for 
the  12  hour  average  wind  speeds).  This  is  further  illustrated  in  Fig.  9 
where  measured  Vs(r)  spectra  are  compared  to  the  WMK  model  prediction  (taken 
from  Eqs.  2,  4  and  5)  along  with  a  similar  type  of  model  from  Fitzgerald 
and  Ruskin  (1976).  The  error  bars  represent  standard  deviations  (again, 
about  a  factor  of  three)  about  the  equilibrium  value. 

The  deviations  from  average  and  equilibrium  values  are  more  important 
for  operational  usage  and  somewhat  less  important  for  climatological  and 
spectrum  evaluation  usage.  An  example  where  both  usages  are  affected  is  the 
estimation  of  1R  extinction  using  visible  extinctions  (visibility  observations) 
and  graphs  similar  to  Fig  4.  We  have  replotted  the  visible  to  IR  ratio 
in  Fig  10  with  lines  indicating  the  WMK  model  predictions  using  two  dif¬ 
ferent  continental  aerosol  coefficients.  Since  the  continental  aerosol 
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coefficients  are  correlated  to  visibility  (Fig.  6),  the  correct  ratio 
depends  not  only  of  wind  speed  but  also  the  visibility  observation. 

Although  we  have  dealt  with  the  stochastic  properties  of  ensemble 
averages  of  the  aerosols,  variations  about  the  average  are  not  necessarily 
random  but  are  primarily  due  to  changes  in  synoptic  and  mesoscale  weather 
patterns.  In  the  case  of  the  continental  aerosol  component,  this  is 
basically  a  question  of  air-mass -history.  In  the  case  of  the  sea  salt 
aerosol,  it  is  a  question  of  changing  surface  generation  rates  (wind  speed) 
and  the  production,  removal  and  mixing  mechanisms  in  the  marine  atmospheric 
boundary  layer.  Since  the  surface  generated  aerosols  are  quickly  mixed 
vertically  to  fill  the  boundary  layer,  rapid  changes  in  the  boundary  layer 
height,  h,  will  be  reflected  in  changes  in  the  sea  salt  aerosol  density 
and,  therefore,  the  extinction  coefficient.  In  Fig  11  we  can  see  that 
fractional  variations  in  h  are  highly  correlated  with  variations  in  the 
visible  extinction.  The  correlation  with  10.6  urn  extinction  is  considerably 
less  because  the  large  size  aerosols  (which  are  heavier  contributors  to  IR 
extinctions)  reach  equilibrium  more  quickly  after  changes  ir.  surface 
conditions. 

The  findings  and  recommendations  of  this  study  are  summarized  below: 

1)  The  WMK  continental  coefficient  should  be  changed  from  B  =  1.7 
to  B  =  0.24 

2)  For  climatological  purposes  (Fig  3),  the  model  predictions  of 
average  extinction  are  very  good  (12  hour  average  wind  speeds 
preferred) . 

3)  For  operational  purposes  (Fig  5),  the  model  predictions  of 
observed  extinction  systematically  overestimate  the  extinction 
under  high  visibility  conditions.  The  random  error  is  typically 
a  factor  of  3. 

4)  Calculation  of  IR  extinction  from  visibility  estimates  (Fig  10) 

is  subject  to  systematic  error  due  to  variations  in  the  continental 
background  (Fig  6). 

5)  Changes  in  visible  extinction  over  four  hour  periods  can  be  highly 
correlated  with  changes  in  the  mixed  layer  height. 
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Meteorological  and  Aerosol 
Data  from  JASIN.  The  Time 
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APPENDIX  B 


Meteorological  and  Aerosol  Extinction 
Data  from  CEWCOM-78.  The  Time  is  PDT. 
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